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Abstract 

We derive expressions for the neutrino mixing parameters that result from com¬ 
plex perturbations on (1) the Major ana neutrino mass matrix (in the basis of charged 
lepton mass eigenstates) and on (2) the charged lepton mass matrix, for arbitrary 
initial (unperturbed) mixing matrices. In the first case, we find that the phases of 
the elements of the perturbation matrix, and the initial values of the Dirac and Ma- 
jorana phases, strongly impact the leading-order corrections to the neutrino mixing 
parameters and phases. For experimentally compatible scenarios wherein the initial 
neutrino mass matrix has pi — t symmetry, we find that the Dirac phase can take any 
value under small perturbations. Similarly, in the second case, perturbations to the 
charged lepton mass matrix can generate large corrections to the mixing angles and 
phases of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix. As an illustration 
of our generalized procedure, we apply it to a situation in which nonstandard scalar 
and nonstandard vector interactions simultaneously affect neutrino oscillations. 
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1 Introduction 


After decades of neutrino oscillation experiments, the mixing pattern in the lepton sector 
has been well established [1]. There are one small and two large mixing angles, and two 
mass-squared differences that differ by a factor of 30 in the neutrino sector. Numerous 
neutrino mixing scenarios have been proposed in the literature to explain such a nontrivial 
mixing pattern; for a recent review, see Ref. [2]. The most attractive scenarios are those 
with mixing patterns motivated by simple symmetries, such as tri-bimaximal mixing [3], 
bimaximal mixing |1], and golden ratio mixing [S]. All three mixing scenarios have ^13 = 0, 
623 = 45°, and are a subset of the more general j-i — t symmetry |6l [7]. However, recent 
measurements of ^13 by short-baseline reactor experiments Daya Bay |8], RENO |9], Double 
Chooz [in], and long-baseline accelerator experiments T2K [TT], MINOS [12] strongly disfavor 
013 = 0. Therefore, models with simple symmetries need additional features to explain the 
observed neutrino mixing pattern. 

A modified approach to explain the data is to treat the simple mixing scenarios as the 
underlying model and add perturbations to accommodate the discrepancy between theoret¬ 
ical predictions and experimental data. In Ref. [7] , we took fi — t symmetry (in the charged 
lepton basis), as the underlying model, and added real perturbations to the Majorana neu¬ 
trino mass matrices to explain the data. We found that small perturbations can cause large 
corrections to 0i2, and the experimental data can be explained by most /i — r symmetric 
mixing scenarios with perturbations of similar magnitude. 

After the discovery that the mixing angle 0i3 is relatively large, many scenarios without 
fi — T symmetry have been proposed [13] • Motivated by this, and by the fact that real 
perturbations have no effect on the Dirac and Majorana phases, in this paper we consider 
arbitrary initial mixing for the underlying model and generalize the perturbation results to 
the complex space. Under the assumption that the charged lepton mass matrix is diagonal, 
we derive analytic formulas for the leading-order (LO) corrections to the three mixing angles 
and the Dirac and Majorana phases. We find that the phases of the elements of the pertur¬ 
bation matrix, and the initial values of the Dirac and Majorana phases, strongly impact the 
LO corrections to the neutrino mixing parameters. We also perform a numerical study of 
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complex perturbations on initial nentrino mass matrices with /r — r symmetry. We find that 
the Dirac phase can take any valne under small perturbations for experimentally compatible 
scenarios. 

Since the mixings in both the charged lepton sector and nentrino sector contribute to the 
observed Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, we explore the case in which 
the charged lepton mass matrix is not diagonal, and consider small complex perturbations 
to the charged lepton mass matrix as well. Since the initial mixing matrix in the charged 
lepton sector is nnconstrained, small perturbations in the charged lepton sector could have 
large effects on the 1-2 mixing in the charged lepton sector, which lead to large changes in 
the mixing angles and phases in the PMNS matrix. 

In addition, as an application of our generalized pertnrbation procednre, we study neu¬ 
trino oscillations with matter effects from both nonstandard scalar and nonstandard vector 
interactions. Nonstandard scalar interactions add small perturbations to the nentrino mass 
matrix, which yield corrections to the vacnnm mixing angles and mass-sqnared differences. 
By using our formalism, we demonstrate how expressions for nentrino oscillation probabili¬ 
ties that simultaneously depend on nonstandard scalar and nonstandard vector interactions, 
can be obtained. 

This paper is organized as follows. In Section 2, we work in the diagonal charged lepton 
basis and derive analytic formulas for the mixing parameters that result from real and com¬ 
plex perturbations. In Section 3, we perform a numerical analysis of complex pertnrbations 
to nentrino mass matrices with /i — r symmetry. In Section 4, we calculate corrections to 
the three mixing angles and phases in the PMNS matrix from pertnrbations in the charged 
lepton sector. In Section 5, we apply our perturbation resnlts to study neutrino oscilla¬ 
tions with both nonstandard scalar and nonstandard vector interactions. We snmmarize onr 
results in Section 6. 


2 Perturbations on the neutrino mass matrix 

Onr goal in this section is to obtain the LO corrections to all the physical parameters under 
small pertnrbations on the initial Majorana nentrino mass matrix assnming the charged 
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lepton mass matrix to be diagonal. The final (resnltant) mass matrix can be written as the 
snm of an initial matrix Mq and a pertnrbation matrix E, i.e., 


M = Mq + E 


u;mouI + 


h 

^12 


ei2 

^22 

€23 

^^ei3 

€23 

€33 J 


( 1 ) 


where Mq = diag(m5, mg), and Uq is the initial mixing matrix. 
The final mass matrix can also be written as 


M = U*MU\ (2) 

where U and M have the same form as Uq and Mq. From nentrino oscillation experiments 
we know that mi and m 2 are nearly degenerate, so here we assnme |eij|, |5m2i| 


2.1 Real case 

For simplicity, we consider the real case first. The mixing matrix Uq can be written as 

f/0=i?M3^?2, (3) 

where is the rotation matrix in the i — j plane with a rotation angle 6 ^,. Then Eq. ([T]) 


can be rewritten as 


M = m?I + i?M3^?2 


0 0 0 

0 6 m 2 i 0 




(^?2)^(^?3)^(^23)^ + ^ 



V° 

0 

<^”^31/ 




^0 

0 

0 ^ 


m?I + 3R°33R?3 

Ri2 

0 

5ml^ 

0 

+ E' 



1° 

0 






(4) 
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where E’ = E6m^- = m° — m°, and I is the 3x3 identity matrix. We 

employ the following notation: 

ei = eii , 62 = ei2C23 — ^13523 , €3 = €12523 + ei3C23 , 

€4 = €22(€23)^ + €33(523)^ - €235^23 ! 

£5 = e 23 C 2 x 23 + 2^*'22 “ €33)52x23 ) 

€6 = €22(523)^ + €33(c 23)^ + €2352x23 ) (5) 

where c°-, 5 °-, c^xj^ and 52 xij denote cos 6 *°-, sin 6 *°-, cos(26*°.) and sin(26*°.), respectively. Then 
= {E')ij can be written explicitly as 

^11 = ^1(^13) +^6(5?3) —€352x13, ^12 = ^2C?3 — €555*3 , 

^13 = ‘^3C2x13 + ~ <^6)52x13 > ^22 = <^4 , 

€23 = € 25°3 + € 505*3 , €33 = € 1 ( 55 * 3 )^ + € 6 ( 05 * 3 )^ + € 352 x 13 • ( 6 ) 


In order to obtain the hnal mixing matrix that diagonalizes M, we use a procedure that is 
similar to that in Ref. [H] . We hrst put zeros in the 2-3 and 1-3 entries of the matrix in the 
square bracket of Eq. (|1]) by using rotations R 23 {S 2 ^) and i?i3(553), respectively. To LO in 
0{\eij\/\5rn!l]\), we have 


O 23 ~ 


-23 


(5mo 


31 


0l3 ~ 


-13 


(5mq 


31 


and the LO correction to rUg is 


(7) 


5m3 = €33 . 


( 8 ) 


Note that since \eij\ <C after the two rotations in the 2-3 and 1-3 planes, the matrix 

in the square bracket of Eq. (|1]) becomes block diagonal and the 1-2 submatrix remains 
unchanged to leading order. Hence we can rewrite Eq. (jl]) as 


M = m\\ + V 


M' 0 


0 

where R = R^23Ri3R23{^23)RM3)Ri2, and 


R^ + 0(|€.,f/<5m°3), 


M' = 


^n('^12)^ + ^22('5 i2)^ ^12'®2x12 ^12*^2x12+ 2^^11 ^22)'52x12 


(9) 


f’ r' 
6126 


'2x12 + 2^^11 ^22)'®2xl2 ^il('®12)^ + ^22('-12)^ + ^i2'52xl2 + ^^ 21 ^ 


( 10 ) 
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which can be diagonalized by the rotation Ri 2 {^') with 



,0 

2x12 


+ 


( 11 ) 


The corrections to mi and m 2 can be written as 



( 12 ) 


where A = + 4 (e; 2)2 + “ f'n)^ + 2 i 5 ™ 2 l [ 2 t'i 2 S 2 xi 2 + (42 “ 4 l) 4 xi 2 ]. and the plus 

(minus) sign is for i = 1 (2). The final mass matrix is diagonalized by the following mixing 
matrix 


U = R«AR,3(S',,)Ra(Su)Iii2Ri2(i') ■ 


( 13 ) 


By comparing it to the standard parametrization, we find the LO corrections to the three 
mixing angles to be 



( 14 ) 


where we have ignored the next-to-leading-order correction to 6 ^ 12 , which is 0{\eij\/\5m\i\). 


For 05*3 = 0 and 6 ^^ = 7 r/ 4 , it is easy to verify that the corrections in Eq. (ITT)) yield the results 
of Ref. [7] for the LO correctionsj^ which were obtained using degenerate perturbation theory. 
As noted in Ref. [7], the near degeneracy of mi and m 2 (|^^ 2 il implies that 6612 

can be large for small perturbations (|ejj| -C l^mg^l). 

2.2 Complex case 

For the complex case, the most general form for Uq is 


Uo = R23{0^23)Uis{9^i„6^)Ri2{9^i,)P{4,4), 


(15) 


^Also, for the next-to-leading-order correction to 012, we obtain Eq. (14) of Ref. [7], except that 6 m 2 i in 
the denominator should be replaced by Sm^i = (m^ -t Sm 2 ) — (jn\ + Smi). 
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where 



(1 

0 

0 ^ 



/ „o 

013 

0 6- 

-i5°„0 \ 
*13 

fi23(«S3) = 

0 

gO 

023 

„o 

*23 

, £'i3(e;3.'5“) = 


0 

1 

0 


(0 - 

-sO 

*23 

^23 J 



\ ^ *13 

0 

C?3 ) 


( 4, 

„0 

*12 



^1 0 

0 ^ 


RMi) = 

„0 „0 
*12 O32 

0 


0 g*02/2 

0 



( 0 

0 

1) 


\0 0 

coo 

to 



Because of the nonzero Majorana phases, in general, the mixing matrix would not remain 
unchanged if we subtract the identity matrix multiplied by a constant from the mass matrix. 
Hence we use a slightly different procedure to obtain the LO corrections for the complex 
case. We rewrite the hnal mass matrix as 


M = u;mouI + E = U* 


Mq + E 


uL 


where E = Uq EUq can be explicitly written as 


/ 


E = 


a 














with 


(17) 


(18) 


a = ei{si 2 )‘^ + [ei(c'( 3 )^ — ](ci 2 )^ + — e2C?3)s2xi2 ^ 

b = e2C?3C2xi2 + [ei(c'j’3)^ - 64 + e6(s?3)^e^*^ ]c?2Si2 “ ['^3S2x12'Si3Ci 3 + ^sC^x 12'Sl3]6*^ ) 
C = 64 ( 052 )^ + [ei(c53)^ - 6352x13^*'^ + 60 ( 553 )^ 6 ^*'^ ](’S?2)^ “ (e5-S?3e*'^ - e2c5’3)s2xl2 ) 

d = (eiC° 2 C ?3 - e 2 S° 2 )s? 3 e"*'^° - e 6 C° 2 C? 3 S° 3 e*'^° + e3C?2C^xi3 “ (^ 5013^2 , 
f = (eis52C?3 + 62052)5336 — ^6Si2C%Si^^ + 63s52c5x 13 + 05c53c52 , 


9 — + €382x13^ +'^i(' 5 i 3 )^e 


0 -i5° 


..0 \2-2i5° 


( 19 ) 


Similar to the real case, we apply a unitary matrix Us to N = Mq + E such that there 
are zeros in the 2-3 and 1-3 entries of the matrix UjNUs- Since \eij\ |(5m5i|, to LO in 
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0{\eij\/\5m\^\), Us can be written as 


where 


Us 


1 

0 


0 

1 



<^23 5 


1/ 


(5i3 


mg — 1 ’ 


^23 ~ 



mg — m5e“^*‘^231 ’ 


(20) 


( 21 ) 


with tan^ig = tan [arg(ci) + 0 ^/ 2 ] and tan ())23 

block-diagonalization, the LO correction to m 3 is 


mg+m 
777 . 3 —m 


0 

1 

1 


tan 


arg(/) + 



After 


(5m3 


mg + fife’ 


— m 


0 

3 • 


( 22 ) 


Note that the 1-2 snbmatrix of N remains nnchanged to leading order after the block- 
diagonalization. Using the procednre described in Appendix|^ we diagonalize this snbmatrix 
nsing the nnitary matrix 


where 


Ul 2 {iA) 


cg 0 

V 0 0 V 


(23) 


, ^ \a + mP^\sm{(j)a 

0 = arctan -----— 

I a + mp cos( 0 (j 


(j)b) — + m^l sin( 0 c — </>&) 

(j)b) + |ce *‘^2 777.2 I cos(0c “ 4>b) ’ 


(24) 


f = - arctan ■ 
^ 2 


2 \b\ 


(25) 


|ce *'^2 - 1 _ cos(0c + 0 — 06) — |a + n7°| cos(0a — 0 — 0b) ’ 

with 0a = arg(a + m?), 0b = arg(&) -f 02/2 and 0c = arg(ce *'^2 + 772 °). In addition, we obtain 
the LO corrections to mi and m 2 as 


dull = 
51712 = 


(a + m?)c| + (ce*'^° + m°)s|e 2 *^ - 2 bs^c^U^ 
(a + m?)s|e“^*'^ + (ce*'^° + m2)c| + 2bs^c^e~ 


m 


1 ) 


— m 


0 

2 • 


(26) 
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The final mixing matrix that diagonalizes M and makes the diagonal elements real and 


non-negative can be written as 


where P = diag(e*‘^i/^, and 

ui = — arg (a -I- m°)c| -f (ce *'^2 -|_ 77 i°)s 2 g 2 *</< _ 2bs^c^P' 

0 J 2 = - arg (a -f + mDc^ + 2 bs^c^e~'"^ 

cjs = — arg -|- j . 


(27) 


(28) 


As shown in Appendix [HI the right-mnltiplication of Ui 2 {^,(f>) does not change 6*13 and 
623 - Hence, the LO corrections to 613 and 623 come from the right-mnltiplication of Us- 
Since ^13 and ^23 are snppressed by a factor of lejjl/l^mggl, while ^ and 0 are not, the LO 
corrections to Ou and the Dirac phases come from the right-mnltiplication of Ui 2 {^, 4 >), and 
the LO corrections to the Majorana phases come from both Ui 2 {^, 4>) and P. 

By comparing U to the standard parametrization, we obtain the LO corrections to the 
three mixing angles: 


rn _ M|Ci2 COs(( 5° - ^ - 013) |/|s?2 COs((5° - ^ 23 ) 

<jtfl3 — -1- n -n-1-T 


jmg — 1 


[mg — 


( 29 ) 


\d\42 COs{^ + 4 >i3) |/|c;2COs(^h^ -fe) 

0u23 — -;— 1 \ -f;———; -r 


jmg — 1 


[mg — 1 


(30) 


59 i 2 = arcsin \ sin^( 6*)’2 + 0 “ sin( 26 *j' 2 ) sin( 2 ^) 


•2 02 + 20 Q 


sm 


^ 12 , 


where t?,- denotes tanAk. The LO corrections to the three phases can be written as 

tjll 




A6 = a — 0 , 


( 31 ) 


(32) 


- —2{a + /I) + W2 - Wi, 


(33) 


A03 — -2/1 + . 


(34) 
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where 


tan^?, 2 + 6 ) 

a = — arctan---- tttt, -w 5 

1 — tan6'°2 tan.^ cos(02/2 + 0) 

and 


( 35 ) 


B = arctan tanesin(0°/2 + 0) 

^ tan 0?2 + tan ^ cos(0O/2 + 0) ' ^ ^ 

From Eq. flHT]) . we see that 8612 varies from to +.^ depending on the initial Majorana 
phase 02 and the pertnrbation phase 0 . Since ^ and 0 depend only on the ratios of linear 
combinations of e^-’s and large corrections to 6*12 and the Dirac and Majorana phases 

are possible even for small pertnrbations. However, corrections can be small in special cases, 
e.g., if 02 is close to 180 ° for the inverted hierarchy, 0 approaches 90 ° and ^ is snppressed 
by a factor of |eqj/(m2 + m^), so that the corrections to 612 and the Dirac and Majorana 
phases are also small. 

Note that the corrections in the complex case are strongly dependent on the phases of 
ejj, and the initial valnes of the Dirac and Majorana phases. If we take e^’s to be real, and 
set ( 5 ° = 02 = 03 = 0 in Eqs. (I 29 |) . (I 30 D and (l 3 TD . we recover Eq. (IT^ . 


3 Perturbations to /i — r symmetry 

As an illnstration of onr analytic resnlts, we stndy pertnrbations on initial nentrino mass 
matrices with p — r symmetry. There are fonr classes of mixing with /i — r symmetry [7j: (a) 

^23 = 45°,0?3 = 0 ; (b) = 45°,0?2 = 0 ; (c) 0°3 = 45°,0?2 = 90 °; (d) 6^^ = 45 °, = ± 90 °. 

In Ref. [ 15 ], it was shown that the initial class (a) can be pertnrbed to class (d) for a specihc 
model. Here we reprodnce the resnlts of Ref. [T 5 | by applying onr general pertnrbation 
formnlas. The complex nentrino mass matrix of Ref. [H] can be written (in onr phase 
convention) as 



^1 ± 25 '" 

0 

0 ^ 


( 

25 ' 

5 " 

_5".\ 

M = mo 

0 

5 '" 

-(1 ± 5 '") 

± rno 


5 " 

0 

0 



1 0 

-(1 ± 5 '") 

5 '" y 


V 

- 5 "* 

0 

0 

/ 
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where ttiq is a common mass parameter, 5'", 6 ' are real and |5'|, |(5"| We treat the 

hrst term on the right-hand side of Eq. fn7|) as the initial mass matrix and the second term 
as the perturbation. The initial mass matrix has class (a) /i — r symmetry. In the standard 
parametrization, we have = j, 6^2 = ^13 = 02 = 0 03 = The three initial masses 

are = ^ 0(1 -|- 26'"), and = rriQ. In this case, Eq. fll9p is greatly simplihed: 


a = 2mQ6 ', b = \/2moRe {5") , c = f = g = 0 , d = i\/2moIm (5") . (38) 


From Eqs. fl2lll . fl2Tl) and fl25ll . we hnd 

Im (5") 


^23 — 0 , 


Jl3 


V26'‘ 


= 0 , 


^ 1 ^ x/2Re (5") 

^ 2 - 6 ' 


Then the hnal mixing matrix can be written as 




U = R 23 [-)P( 0 , 7 r)R„[Su)Ru(() 

= R2^(J)Ua(^a,^)RMi)P(0,’')■ 


( 39 ) 


( 40 ) 


Hence, the hnal mixing angles and the Dirac phase are 


6*23 — 


TT 

4 ’ 


6*12 = - arctan 


V 2 Re 


6*13 — 


Im ((5" 

V25" 


6 = 


TT 


(41) 


as in Ref. na. Note that the initial class (a) is perturbed to class (d), and that the large 
change of the Dirac phase S coincides with the deviation of ^13 from 0. 

The general form of the neutrino mass matrix with class (d) fi — t symmetry and its 
associated generalized CP symmetry has been recognized in Ref. [16], and deviations from 
it were discussed in Ref. na. It has been shown in Ref. [IB] that the general forms of the 
neutrino mass matrices with class (a) and (d) /i — r symmetry are (in our phase convention) 


(42) 


respectively. Here x, y, z, w, r, s are complex and u, v are real. Hence, any perturbation 
matrix of the form 



X 

y 

-y 


u 

* 

1 

Ma = 

y 

z 

—w 

, and Md = 

r 

S —V 


\-y 

—w 

^ J 


[-r* 

—V S* i 


/ 


E = 


-12 


-f 2 ilm{y) 


Re(eii) - Hm(a;) ei 2 

ei2 622 Re(e23) + Hm(t(; 

-e ^2 + 2 ilm{y) Re(e 23 ) + Hm(u;) e ^2 “ 2 i\m{z) j 


(43) 
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Table 1: Best-fit values and 2 a ranges of the oscillation parameters [IB], with 


and Am^ = — (mf + ml)/ 2 . 


Parameter 


9i3(°) 

023(°) 


Am^ (10 ^eV^) 

Normal hierarchy 

33.7l|} 

8.801°;^? 

41.416-6 

7 KA+0A& 

' •'-'^-0.39 

2 43+0-12 

Inverted hierarchy 

33.7l|} 

8.9ll°i° 

42.4I|6 

7 KA+0A& 

' •"-'^-0.39 

0 00+0.12 
+ -'l°-0.13 


perturbs the initial mass matrix with class (a) p — r symmetry to class (d) /i — r symmetry. 

We now perform a numerical search to find perturbations that fit the experimental data 
for initial neutrino mass matrices with /i — r symmetry. We select class (d) and scan 6^2 
and 9^^ over the range [0,90°]. Since the initial mass matrices of classes (a), (b) and (c) 
do not depend on 5°, the perturbation results of class (d) will cover the other classes, e.g., 
the perturbation results for bimaximal mixing would be the same as that of class (d) with 

= 0 and 6^2 = 45°. Since we work in the basis in which the charged lepton mass matrix is 
diagonal, the mixing matrix in the neutrino sector is the same as the observed PMNS matrix. 
We also choose mi = 0 for the normal hierarchy (or m 3 = 0 for the inverted hierarchy), so 
the best-fit values from the global fit in Table [T] define the other two final masses and the 
three hnal mixing angles. 

We characterize the size of the perturbation as the root-mean-square (RMS) value of the 
perturbations. 


Crms — 


Tr[EtE] 



(44) 


where i and j sum over neutrino flavors. eRMS is determined by the three initial masses, two 
initial Majorana phases, two final Majorana phases and one final Dirac phase. 

The initial Dirac phase in class (d) is fixed to be ±90°. To evaluate the change in the 
Dirac phase due to the perturbations, we £x 6^2 = 45° and scan over 6 and 6 ^^ to hnd the 
minimum RMS value of the perturbation, that results in the best-fit parameters. The 

results for 6 ^ = 90° are shown in Fig. [1] The results for 6 ^ = —90° (or 270°) are symmetric 
to those of (5° = 90° with 6 —)■ 360° — 6 . From Fig. [H we see that for 6 ^^ < 20° it is possible 
for the final Dirac phase to have any value under small perturbations (e]^^g < 10 meV), i.e.. 
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Figure 1: Iso-eg';J5g contours in the ( 6 *° 3 ,( 5 ) plane that give the best-ht parameters for /i — r 
symmetry with ^23 = ^12 = 45°, <5° = 90°. The left panel is for the normal hierarchy with 
mi = 0 and the right panel is for the inverted hierarchy with m 3 = 0 . 

the correction to the Dirac phase can be large for small perturbations. 

4 Perturbations in the charged lepton sector 

In the basis in which the charged lepton mass matrix is not diagonal, the observed PMNS 
mixing matrix is 

Upmns = U}U,, (45) 

where Ui and are the mixing matrices in the charged lepton and neutrino sectors, respec¬ 
tively. For an arbitrary charged lepton mass matrix Mi, we have 

{Mi)^Mi = UM\Uiy, (46) 

where Mi = diag (mg, m^, mr). 

Snppose the charged lepton mass matrix is also the result of small perturbations to an 
initial mass matrix, i.e.. Mi = -|- Ei, where {Ei)ij = eh and |eh| mr- If the initial 
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mixing matrix in the charged lepton sector is 17°, i.e., 




( 47 ) 


then to LO, we get 




= U? 






(48) 


where = {U^Y [{Mj^YEi + 17°. Note that since 17° is unconstrained, the size of 


each element of the iV^ matrix could be of order mT-|eh|. 


If (M°)^ + is diagonalized by a unitary matrix 17], i.e., 

Ni An \ 




M° +N^ = 




(NYr Hr + ^22 

I f ]\ yI _i_ atI 


\ (H) 


= ul 


Hr Hf+N^^j 


rUg 0 0 

0 ml 0 
0 0 ml 


\ 


(vW , 


(49) 


then the PMNS matrix can be written as 


UpMNS = {U?Ul)‘'Ut = (Ul)‘'U„, 


(50) 


where Uq = (17°)ll7° has the most general form of Eq. ([T5]). Since ~ m,-|eh|, the 2-3 
and 1-3 mixing angles in 17] are suppressed by a factor of Nljml ~ |eh|/m,-. To LO in 
C>(|e] |/mT-), 17] can be parametrized as 




0 

1 




^23^ *‘^23 


0 —5]3e*‘^23 


/ 


1 

0 


^ [ cos6\ 


0 5\^e *'^13 


y—0 


0 

1 


12 


sm 


sin 6 *] 2 e*'^i 2 cos 0^2 
0 0 


^«2e-*<^i2 o\ 
0 


1 

(51) 


where 


Oi3 ~ 


\Ni 


131 


lA'il 


mt 


->23 


mt 


'12 


1 

- arctan 


2\Ni 


121 


ml + N', - N‘, • 


(62) 
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and (pij ~ — a.TgN-j. 

If 0^2 is also very small, the LO corrections to the three mixing angles in the PMNS 
matrix are 


S0i3 = - 0 I 2 S 23 COs(5'^ - ^ 12 ) - ^13^23 COs(5‘^ - 013) , 


l 


SO23 = -6^3 cos 023 - <^13^23^13 COs(5° - 0i3) + 0i2cL^13 COs((5° - 0i2) , 

SO12 = “o” (^13'®23 013 ~ ^12*^23 0i2) • 


(53) 


n3 


However, in general, since Nlj ~ mT-|eh|, and if |eh| ~ mf^/rrir = 6 MeV, 9[2 conld be very 
large, which will give large corrections to the mixing angles in the PMNS matrix. Thns, the 
sitnation in the charged sector is similar to that in the nentrino sector: the near degeneracy 
of rrie and 771 ^ (on the scale of m,-) can lead to large corrections in 1-2 space. 

For large 6 [ 2 , the analytical expressions for the corrections to the mixing angles in the 
PMNS matrix are cnmbersome. Here, as an illnstration, we consider the very simple scenario 
in which 


/ 


Ul = 


cos ^12 


sin 6 ' 


12 


sin 0(2 
cos 6*12 


V 




(54) 


Then from Eq. 


0 0 

, the hnal mixing angles in the PMNS matrix are given by 


C13C23 — C?3C23 , 


(55) 


4 = ( 4 )^( 4 )^ + ( 4 )^( 4 )^( 4 )^ - 244442^2 cos5°, 


(56) 


'^13'®12 


= [( 


pO 0 


444)^ + (4)^(4)^(4)^(4)^ 


+ 2Si2Si2S?3S23(Ci24'^12 51244)*^°®'^°] ! 


(57) 


where c \2 denotes cos 0 ( 2 , and s [2 denotes sin 6 *( 2 - As an example, if 9{2 is the Cabibbo 
angle and the initial PMNS matrix has bimaximal symmetry (^12 = 45°, = 0), then the 

resnlting ^12 and 9 i 3 are consistent with the observed valnes to within 2 a. 

There are eight parameters in Eqs. ([55]), fl56|l and fl57|) . We nse the best-£t valnes in 
Tabled] for the normal hierarchy to £x 6 ^ 12 , 6*13 and 023- Then for given valnes of ^13 and 5°, 
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Figure 2: Dependence of 6^23 on 6 ^^ for small perturbations in the charged lepton sector when 
Ug is given by Eq. fIMl) . and the three mixing angles in the PMNS matrix are hxed by the 
best-fit values of the global fit in Table [H for the normal hierarchy. 

the other three unknown parameters 6 * 23 , ^12 and 9^2 are determined by the three equations. 
First, we obtain 6*23 from Eq. (I55|) for a given 6 * 33 . Note that the constraints on ^23 and ^33 
are symmetric for fixed ^13 and 623 , which can be seen from Fig. [2l Then we scan 6[2 from 
[—90°, 90°] to hnd solutions to Eq. fl56l) for a given 5°. For each solution of 6 ^ 32 , we obtain 
6^2 from Eq. fl5711 by scanning 6^2 from [—90°, 90°]. Note that we only scan the hrst and 
fourth quadrants of 0^2 [ 6 * 12 ]! because Eq. fl56|) [Eq. fl57|l ] is only sensitive to the relative sign 
between the cosine and sine of 6*^3 [ 6 '? 2 ]- Once we obtain 9 ^ 3 , 9\2 and 6*32 for given values of 
6^33 and 5°, the resulting PMNS matrix is completely determined (except for the diagonal 
Majorana phase matrix) from Eqs. (I50D . (I5T|) and ([15]). By comparing the PMNS matrix 
with the standard parametrization, the resulting Dirac phase 6 is also obtained for given 
values of 6*33 and 5°. The dependence of 6 *^ 2 , 6*12 and 5 on 5° for different values of 6^33 is 
shown in Fig. [3l From Figs. [2] and [3] we see that the initial mixing angles and the initial 
Dirac phase can be very different from their observed values in the PMNS matrix due to 
small perturbations in the charged lepton sector. 

Generally, perturbations in both the charged lepton and neutrino sectors will be present. 


16 




In this case, one must first use the procedure described in this section to hnd the corrections 
to the initial PMNS matrix from perturbations in the charged lepton sector alone, then use 
the new PMNS matrix to rotate to the basis in which the final charged lepton mass matrix is 
diagonal, and ultimately use the procedure described in Section 2 to find the final corrections 
to the parameters in the PMNS matrix from perturbations in the neutrino sector. 


5 Neutrino oscillations with nonstandard interactions 


We now apply our generalized perturbation procedure to a phenomenological study of neu¬ 
trino oscillations that are affected by nonstandard scalar and nonstandard vector interactions 
simultaneously. 

As Vf. propagate in matter, they scatter on electrons via the V-A interaction. This is 
described by the MSW potential |T9], which is added to the vacuum oscillation Hamiltonian: 


E 


‘IE,, 


-U 


mt 


0 0 




V 


0 0 

0 0 ml J 


U^ + 


V 


0 0 

0 0 0 
0 0 0 




(58) 




where Gp is the Fermi constant, W is the electron number density in the medium, and U 
and TTij are the mixing matrix and eigenmasses in vacuum, respectively. 

New physics beyond the Standard Model can be probed by studying model-independent 
nonstandard interactions in neutrino oscillation experiments; for a recent review see Ref. |20j . 
Most studies of nonstandard interactions are focused on the vector interaction, which can 
be described by effective four-fermion operators of the form 


[/Fp( 1 ± 7*^)/] + h.c., (59) 

where / = m, d, e is a charged fermion held, and e)(^ are dimensionless parameters that denote 
the strength of the deviation from the standard interactions. Similar to the MSW term, the 
matter effect due to the nonstandard vector interaction modihes the oscillation Hamiltonian 
by additional potential terms, \/2 GfN 

In addition, consider nonstandard scalar interactions, which may arise from a Lagrangian 
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of the form 

■Cs = , 


(60) 


where 0 is a new scalar held, and A“^ and A/ are dimensionless coupling constants for neu¬ 
trinos and charged fermions, respectively. In a mean held approximation, the nonstandard 
scalar interaction will shift elements of the neutrino mass matrix by 

A“^ 






-\jNf , 


( 61 ) 


where is the mass of the scalar held, Nf is the number density of the charged fermion /, 
which is assumed to be nonrelativistic. 

Tests of the inverse square law of the gravitational force put stringent m^-dependent limits 
on the coupling of a new scalar held to the nucleon held [22]. For in the range, 10“® eV to 
10“^° eV, the current experimental upper limit of Aat varies from 10“^^ to 10“^^ |23|. Since 


— 0.46 meV 


A. 


10 


-4 


A 


N 


10 


-21 


N 


f 


Na/ cm^ 


lO-^* eV 

7710 


( 62 ) 


and Nf ~ liV^/cm^ ~ 10^° eV^ on Earth and Nf ~ lOOA^/i/cm^ ~ 10^^ eV^ in the solar core 
where most solar neutrinos are produced, in these environments, a A^, of order 10“^ gives a 
mass shift of order 1 meV for = 10“® eV. Such values are possible for much smaller 
values of when rrifj, < 10“® eV. 

In the presence of both nonstandard scalar and nonstandard vector interactions, the 
ehective Hamiltonian for neutrino oscillations can be written as 



(i 

0 

0^ 


(e^ 


^er 


0 

0 

0 

+ y/^G pNf 


MM 

T 


lo 

0 

V 


1 V* 

\ ^er 

M-r 

(^rr) 


(63) 


where the ehective mass matrix has the form 



^ mi 

0 

0 ^ 


Un 

ei2 


Meff = U* 

0 

m2 

0 

U^ + 

^12 

€22 

€23 



0 

m^j 



€23 

€33 y 


We apply our generalized perturbation procedure to the study of both nonstandard scalar 
and nonstandard vector interactions. By incorporating the the corrections to the vacuum 
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oscillation parameters (arising from the scalar interaction) into the nonstandard vector in¬ 
teraction formulas, we immediately obtain new formulas for oscillation probabilities with 
both nonstandard scalar and nonstandard vector interactions. Taking the oscillations of 
in long baseline experiments as an example, the result for the survival probability is [ 2 ^ 


p ~ 1 


’2x23 


sm 


2 

4E 




COS0)(^S2x23 


4x23i^GFNeL) si 


AmliL 


sm ■ 


+ Ac: 


2\/2GFNeE . 2 AruoiL 


2x23 


Amil 


+ (i^: 


V I 


-|)'S2x23'^2x23 


2E 

^/2GFNeL . Aml^L ^2^/2GFNeE 


sm 


AE 


sm ■ 


- 2 - 


sm 


2 AmiiL 


, (65) 


2 2E Amg;^ AE 

where = m\ — ml, S 2 xij = sm29ij, C 2 xij = cos29ij, and (j)^^ = arge^^. After replacing 

Am|g —)■ Am|g -|- 2{mz5mz — mi6mi) and 6*23 —)■ 6*23 + 6923, where the shifts in and 6^23 
can be easily obtained from our perturbation results in Section 2, the new formula for both 
nonstandard scalar and nonstandard vector interactions is as follows: 


— J- *2x23 


sm 


,2 Am|^L 
AE 


—26 923 sin 4023 sin 


2 Am^L [m 36 m 3 — mi6mi)L 


Am\^L 


AE 


2E 


S2x23Sin- 
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~ F/irl cos §2x23 


t^T- 


+ (141 - I4|)^L23C2 


slx23i42GFN,L) 


Aml^L 


sm 


2E 


+ Ac: 


2\/2GFNf,E . 2 Am\,L 


2^23 Amii 


sm 


AE 


x23 


y/2GFNeL . ArUgiL 2\/2GFNeE . 2 Am\,L 


sm 


2E 


- 2 


Amli 


sm 


AE 


( 66 ) 


We see that cancellations between the nonstandard scalar and vector terms are possible, a 
study of which is beyond the scope of this paper. 


6 Summary 

We introduced a generalized procedure to study complex perturbations on Majorana neutrino 
mass matrices. In the charged lepton basis, we derived analytic formulas for the corrections 
to the three mixing angles, and the Dirac and Majorana phases for arbitrary initial mixing. 
Since mi and m 2 are nearly degenerate, the corrections to 0 i 2 and the Dirac and Majorana 
phases could be very large. We performed a numerical analysis on the mass matrices with 
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/i — r symmetry to illustrate our analytical results, and found that the final Dirac phase can 
take any value under small perturbations. 


We also studied the scenario in which the charged lepton mass matrix is not diagonal, 
and considered perturbations on the charged lepton mass matrix. We found that small 
perturbations in the charged lepton sector give small mixing in the 1-3 and 2-3 sectors, but 
the mixing in the 1-2 sector could be potentially large due to the near degeneracy of rUg and 
rrin (on the scale of rrir), which could lead to large corrections to all three mixing angles in 
the PMNS matrix. 

In addition, we showed that using our generalized perturbation procedure, it is straight¬ 
forward to study neutrino oscillations with both nonstandard scalar and nonstandard vector 
interactions. 

Acknowledgements. We are grateful to R. Sawyer for useful discussions. This research 
was supported by the DOE under Grant No. de-sc0010504 and by the Kavli Institute for 
Theoretical Physics under NSF Grant No. PHYll-25915. 


A Diagonalization of a 2 x 2 complex symmetric matrix 


We show how to diagonalize a 2 x 2 complex symmetric matrix 



(67) 


SO that 



( 68 ) 


where mi, m 2 are non-negative real numbers, and U is an unitary matrix. 
First, we diagonalize M with a unitary matrix V of the form 



( 69 ) 
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where and denote cos.^ and sin.^, respectively: 

yrpMV I cs|e^*'^ — 2bs^c^e^^ (ae“*'^ — ce*'^)s^c^ + 6(c| — s|) 

y(ae“*‘^ — ce*'^)sgcg + 6(c| — s|) as|e“^*'^ + cc| + 26s^c^e“*'^ 

The diagonalization condition is 

(ae”**^ — ce*‘^)s5C^ + 6(c| — s|) = 0 , 

which implies the phase 0 is 

0 = arctan ■ 

and the mixing angle is 


a 

1 sin(0a - 0fe) - 

c 

sin(0c - 0b) 

a 

cos(0a - 0b) + 

c 

cos(0c - 0b) ’ 


1 


2\b\ 


^ — _3jrct3jH_ 

2 |c| cos(0c + 0 - 0fe) - |a| cos(0a - (f) - (ph) 

Also, the two eigenvalues can be written as 

mi = |ac| + cs|e^*'^ — 26sgCge*‘^| , 
m2 = + cc| + 26s^c^e“*'^|. 

The hnal unitary matrix that diagonalizes M is 

ei^i/2 Q 


U = V 


0 




where 


cji = — arg (ac| + — 26s^c^e*'^) , 

U 2 = — arg + cc| + 2 bs^c^e~'^'^) . 


( 70 ) 


( 71 ) 


( 72 ) 


( 73 ) 


( 74 ) 
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B Right-multiplication 

We calculate the change of mixing parameters when a general initial mixing matrix Uq (see 
Eq. [T5|) is multiplied by the following unitary matrix from the right: 


Ul2{iA) = 


s^e 

0 0 1 


( 77 ) 
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Multiplying Uq on the right by U 12 yields 
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A 
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^ '-13'-23 


. (78) 


where 


and 


(/>2 +2(/> 


<^12 = - s\2 -HR 2 




>§12 = <Sl2C^ + c?2-s^e* "2 
are complex. Comparing U to the standard parametrization, we hnd that 


^23 — ^23 > ^13 — ^13 ) 


and 


612 = arcsin(|^i 2 |) = arcsin \ sin^( 6'°2 + 0 ~ sin( 26 *° 2 ) sin( 2 ,^) sin 


2 02 + 20 


(79) 


(80) 


(81) 


(82) 


Note that after the right-multiplication the phases of the resulting mixing matrix are not 
in the standard form. Dehning 

tan 9^2 tan,^sin( 02/2 -f- 0 ) 


a = arg(Ci 2 ) = — arctan 


0 = arg(Si 2 ) = arctan 


1 — tan 9^2 tan,^ cos( 02/2 -|- 0 ) ’ 

tan,^ sin( 02/2 -f- 0 ) 
tan6'°2 + tan,^ cos(02/2 -f- 0) ’ 


we can write U as 
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U = 


„0 1/^ \„ia 
Ci 3 |Gi 2 |e 

C23l>§i2|e*^ - S23^?3l<^12|e*^^°+“^ 
S 23 l>§i 2 |e*^ - C 23 S? 3 |C'i 2 |e*(^°+") 
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^ 1^13*23 
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Q ^yj ^yj 

^ '-13'-23 


(85) 
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On removing the unphysical phases (j)e = a and (j)^ = (j)r = [3 from the rows, the phases 
in the second and third columns match the form of the standard parametrization, with the 
Majorana phases shifted by 


A02 — —2{a + (3), 


( 86 ) 


A .).3 = - 2 /?, 


( 87 ) 


and the Dirac phase shifted by 


AS = a — f3. 


( 88 ) 
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